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<^> ' We compute the supersymmetric (SUSY) contributions to v (z^)-nucleus deep in 

(N ' , ' . , . , 

elastic scattering in the Minimal Supersymmetric Standard Model (MSSM). We con- 

sider the ratio of neutral current to charged current cross sections, R u and i?p, and 

m , 

£S) , compare with the deviations of these quantities from the Standard Model predictions 

implied by the recent NuTeV measurement. After performing a model-independent 

l ; . ^ - „ r 

the NuTeV anomaly. We discuss one scenario in which a right-sign effect arises, 



o 
m 
O 



X 



and show that it is ruled out by other precision data. We also study for R parity- 
violating (RPV) contributions. Although RPV effects could, in principle, reproduce 



Ph. the NuTeV anomaly, such a possibility is also ruled out by other precision electroweak 

measurements. 
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I. INTRODUCTION 

Neutrino scattering experiments have played a key role in elucidating the structure of the 
Standard Model (SM). Recently, the NuTeV collaboration has performed a precise determi- 
nation of the ratio R u (Rp) of neutral current (NC) and charged current (CC) deep-inelastic 
v n (z^t) -nucleus cross sections jjj, which can be expressed as: 

^(p) = (fff) 2 + r^(flf) 2 , (i) 



where r = /&vN an d (dl r) 2 are effective hadronic couplings (defined below). Comparing 
the SM predictions^, [j| for (gf^j) 2 with the values obtained by the NuTeV Collaboration 
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yields deviations 1 8R v (p) = R c „^ — R^) 

5R U = -0.0033 ± 0.0007, 5R P = -0.0019 ± 0.0016. (2) 

Within the SM, these results may be interpreted as a test of the scale-dependence of the 
sin 2 6 W since the (gl^) 2 depend on the weak mixing angle. While the SM prediction for 
sin 2 6 W at = M z has been confirmed with high precision at LEP and SLC, the predicted 
running of this parameter to lower scales has yet to be studied systematically. The results 
from the NuTeV measurement imply a +3cx deviation at [i ~ 10 GeV, while the current value 
of the cesium weak charge, extracted from atomic parity- violation (PV), implies agreement 
with the predicted SM running at a much lower scale ~\ . Measurements of the P V asymme- 



tries in polarized ee and ep scattering will provide further tests of this running at fi ~ 0.2 
GeV 0Q. 

This interpretation of the NuTeV results has been the subject of some debate. Un- 



accounted for QCD effects, such as charge symmetry-breaking in parton distributions or 
nuclear shadowing?]], have been proposed as possible remedies for the anomaly. Alterna- 



tively, one may consider physics beyond the SM, as reviewed in Ref. |8(. In what follows, 
we focus on one new physics scenario, namely, supersymmetry (SUSY). While a brief dis- 
cussion of SUSY is given in Ref. j^, an extensive, detailed treatment has yet to appear 
in the literature. Because SUSY is one of the most strongly-motivated extensions of the 
SM, undertaking such an analysis is a timely endeavor. The goal of the present study is to 
provide this comprehensive treatment. 

In performing our analysis, we work within the framework of the Minimal Supersymmetric 
Standard Model (MSSM), which remains the standard baseline for considering SUSY effects 
in precision observables. Typically, MSSM studies adopt one or more models of SUSY- 
breaking mediation, thereby vastly simplifying the task of analyzing the MSSM parameter 
space. Here, however, we carry out a model-independent treatment, avoiding the choice of a 
specific mechanism for SUSY-breaking mediation. Since generic features of the superpartner 
spectrum implied by the most widely adopted SUSY-breaking models may not be consistent 
with precision data 9], we wish to determine whether there exist any choices for MSSM 



1 We use the quoted experimental errors on R v and R„, rather than adding the errors on #) 2 in 
quadrature, since the latter are correlated and derived from the experimental cross section ratios. 
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parameters that could account for the NuTeV result, even if such choices lie outside the 
purview of standard SUSY-breaking models. 

We find that it is difficult - if not impossible - to choose MSSM parameters so as to 
improve agreement with the NuTeV result. When R parity is conserved and SUSY effects 
only arise via radiative corrections, the magnitude of their contribution is generally too 
small for the NuTeV anomaly to generate significant constraints. Moreover, for nearly all 
parameter choices, the sign of the SUSY corrections to R u and Rp is opposite to that of 
the NuTeV anomaly. We do find one scenario - not considered in Ref. |8j - under which 
SUSY loops generate a right sign effect with relatively large magnitude, but this scenario is 
presently inconsistent with other precision, electroweak data. 

We also consider possible tree-level contributions from R parity-violating (RPV) inter- 
actions, whose presence would render the lightest supersymmetric particle (LSP) unstable. 
We observe that purely leptonic RPV effects, which arise via the definition of sin 2 9 W , could 
in principle also generate a right sign effect to account for the NuTeV anomaly. However, 
precision data also limits the magnitude of this contribution to be considerably smaller than 
necessary. On the other hand, the NuTeV result does yield new constraints on possible 
semileptonic RPV interactions. 

In short, our qualitative conclusions agree with those of Ref. |8[, though we believe we 
have carried out a more exhaustive analysis. Thus, one would have to look to more exotic 
new physics possibilities, such as a Z' boson coupling only to second generation fermions, 
if SM QCD effects are ultimately unable to explain the NuTeV anomaly. We note that the 
situation here contrasts with that for the PV electron scattering asymmetries, where the 
SM contribution is fortuitously suppressed and where SUSY radiative corrections could, in 
principle, produce observable contributions [lo|. 

The analysis leading to these conclusions is organized in the remainder of the paper as 
follows. In Section II, we discuss general features of MSSM contributions to v ^ (z/ M ) -nucleus 
scattering, including both SUSY loop corrections and RPV effects. Section III gives details 
of the loop computation as well as the scan over the MSSM parameter space. In Section 
IV, we give the RPV analysis, including results of a fit to other precision data. Section 
V contains a summary of our results. Explicit formulae for loop corrections with relevant 
Feynman rules are given in the Appendices. 
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II. v^N (vpN) SCATTERING IN THE MSSM: GENERAL CONSIDERATIONS 

For momentum transfers satisfying \q 2 \ « Mj, the neutrino-quark interactions can 
be represented with sufficient accuracy by an effective four fermion Lagrangian: 

C n NC 

= -^^^7 A (l-75)^E^A[2elP L + 2 e Mg (3) 
= -^ L ^7 A (l-7 5 )^7A(l-75)rf + h.c. , (4) 
where P L = (1 - T5 )/2, P R = (1 + j B )/2 and 

e q L = ll-Q qKv em 2 e w + \ q L (5) 
e R = -Q q K u sm 2 6 W + X R . (6) 

The parameters p^ = p^fi = k v = 1 and \ q L R = at tree-level in the SM. These quantities 
differ from their tree-level values when 0(a) corrections in the SM or MSSM are included or 
when other new physics contributions arise. The precise values of these quantities individ- 
ually are renormalization scheme-dependent. When computing MSSM loop contributions, 



we use the modified dimensional reduction (DR) scheme, in which the spatial dimension of 
momenta are continued into d = 4 — 2e dimensions while the Dirac matrices remain four 



dimensional, as required by SUSY invariance. Quantities renormalized in the DR scheme 
will be indicated by a hat. Note also that for the neutrino reactions of interest here, p„§ , 
p^x, an d k u are universal (independent of quark flavor), while the \ q L R are flavor-dependent. 

The NC to CC cross section ratios R u and R D can be expressed in terms of the above 
parameters via the effective couplings {g C L R ) 2 appearing in Eq. (0) in a straightforward way: 

The SM values for these quantities areQ,[J (gf f = 0.3042 and (gf) 2 = 0.030f while the 
NuTeV results imply (gf) 2 = 0.3005 ± 0.0014 and (gf) 2 = 0.0310 ± 0.0011. 

In what follows, we concentrate on the MSSM contributions to p^v' * 7 , k u , X q L R , R u , and 
Rp. When the R-parity quantum number (— i) 3 ( B - L )+ 2s i s conserved, the MSSM contributes 
to these quantities only via loop effects. The relevant diagrams are shown in Figs. and 
the Appendices. Note that all gauge boson self energy corrections, as well as leptonic vertex 
and external leg corrections, contribute only to the universal parameters p^' CC and k v . For 
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FIG. 1: One loop contributions to the neutrino-quark neutral current amplitude. The blob denotes 
the one loop irreducible diagram or the counter term. The Feynman diagrams for the external leg 
corrections are not shown explicitly. 
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FIG. 2: One loop contributions to the neutrino-quark charge current amplitude. The blob denotes 
the one loop irreducible diagram or the counter term. The Feynman diagrams for the external leg 
corrections are not shown explicitly. 



the NC amplitudes, non-universal box diagrams and quark vertex and external leg correc- 
tions (Syg ,q ) are contained in the X q L R . All CC box graphs as well as hadronic vertex and 
external leg contributions {8y%) appear in p^fi. The CC box graphs also generate ampli- 
tudes involving products of scalar and pseudoscalar currents. However, the corresponding 
0(a) corrections to R v and R p are suppressed by lepton and quark masses, so we neglect 
them here. In terms of these various corrections, the renormalized parameters are 

fi W w(q 2 ) rw(o) 



PvN 



NC 
PvN 



1 + 



1 + 



flzzjq 2 ) 
M 2 z -q 2 ' 



Ml 



+ 8- 



CC ?;u 
VB °VB 



n 



WW 



(0) 



M 2 



+ 8 V V - & l 



VB 



(8) 
(9) 
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cfl z -y(q 2 ) t 5sm 2 w 
8 ? 

Al , R = , (11) 



«, = --^^ + - Ac F A u (q ) (10) 

s q 2 sin 2 9 W 



where liyviq 2 ) are the gauge boson self-energies renormalized in the DR scheme at a scale 
\x = M z ; 5 VB denote vertex, external leg, and box graph corrections entering the muon- 
decay amplitude; and Sy indicates the neutral current neutrino vertex and lepton external 
leg correction. Note that the muon decay corrections enter the semileptonic amplitudes since 
the Fermi constant G M is taken from the muon lifetime. The correction F A ,u arises from the 
charge radius. 2 Superpartner loop contributions to Ilz 7 (g 2 = 0) are zero, so that the 
corresponding contributions to k v are finite at the photon point. The shift in sin 2 9y/ arises 
from its definition in terms of a, G^, and M z 

-2-2 710 

V2M 2 z G,[l-Af(M z )} 



s 2 = l-c 2 = sm 2 6 w {M z ) and|ll| 

a- fr ■ J ^jo) fi zz (M 2 ) , fi ww (o) , ?/t 

Ar = n 77 (0) + 2- t —, ^- + + 6 VB (13) 

with II 77 (g 2 ) = fl 77 (g 2 )/g 2 . Writing Af = Af SM + Af SUSY one has 

^!k§Y = _^_ Af susY (M) 

In Section III, we discuss our computation of the MSSM loop contributions to these parame- 
ters in detail; explicit expression for the various loop amplitudes appears in the Appendices. 

When i?-parity is not conserved, however, new tree-level contributions appear. The latter 
are generated by the B — L violating superpotential 

W RPV = -XijkLiLjE k + X' ijk LiQjD k + -X" jk UiDjD k , (15) 

where Lj and Qi denote lepton and quark SU(2) L doublet superfields, E iy Ui, and Di are 
singlet superfields and the Xij k etc. are a priori unknown couplings. We have neglected 
additional lepton-Higgs mixing term in Eq. (fTK|) for simplicity. In order to avoid unaccept- 
ably large contributions to the proton decay rate, we set the AB ^ couplings X"j k to 



2 The charge radius is equivalent to its anapole moment in the MSSM. 
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FIG. 3: Tree level RPV contributions to muon decay [plot (a)], — q charged current [plot (b)] 
and Ufj, — q neutral current [plots (c) and (d)] amplitudes. 

zero. The purely leptonic terms (Ai 2 fc) contribute to neutrino scattering amplitudes via the 



normalization of CC and NC amplitudes to G M and through the definition of sin 2 9w 



121. 



The remaining semileptonic, AL = ±1 interactions (A- jfc ) give direct contributions to the 
neutrino scattering amplitudes. The latter may be obtained computing the Feynman am- 
plitudes in Fig. El and performing a Fierz reordering. For neutrino-quark scattering, one 
obtains the effective Lagrangian 



T- 



\! 12 \\t 12 

2k1 T u i — 21/r T u i — 



2M k 

I X' I 2 

l A 21fcl 

"r 



a R 



(16) 



where we have taken \q 2 \ « Mj and have retained only the semileptonic terms relevant to 
v^-q scattering. 

In terms of these parameters, the shifts induced in the effective {p^-N parameters 

CC, NC j q 
PvN aIld e L,R are 



$P C UN 

Set 



Am (4) 

S P ^ + A' 21k (d k R ) 

-Kik{d h R ) + -X x A 12k (e k R ) 



5e d R = -A' 2kl (d k L ) + -A,A 12fc (e*) 



5e u L = Se R 



•A,.A ' 



e R ) 



(17) 
(18) 
(19) 

(20) 

(21) 
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where 

IA 



and 
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A^m = p jfc|2 - (22) 



s 2 c 2 

K~~ 2 ^- 0.35 . (23) 

c — s 



The corresponding shifts in R v (p) are 

6R*w = + + r^)A 12k (e k R ) - 2[R% + e d L ]A' 21k (d k R ) (24) 

+2r^e d R A' 2kl (d k L ) 

« -0.25[1 + ^(4) - 2[Rl™ - 0.43]A' 21fc (4) + Lfir^AU^) . 

As we discuss in Section IV, Ai2 k (e R ) and A' 2lk (d R ) are constrained by other precision 
electroweak data, while A' 2kl (d k L ) is relatively unconstrained. In Eq. (J23jl . the coefficients of 
A 21 fe(4) and A 2fel (cf|) are positive, while the coefficient of Ai2fc(e^) is negative. Since the 
Aij k are non- negative, we would require sizable value of Ai2 k (e R ) and rather small values 
of A' 2lk (d R ) and A' 2kl (d k L ) to account for the negative shifts in R v and R v implied by the 
NuTeV analysis. The present constraints on A 12 fc(e^), however, are fairly stringent, ruling 
out sizable values for the semileptonic corrections with fairly high confidence. 



III. SUSY LOOP CONTRIBUTIONS 



Instead of working in a specific SUSY breaking scenario, we perform a model-independent 
analysis by varying . Although such an 

approach is insensitive to the effects of any particular SUSY breaking parameter, it does 
allow us to obtain the size of SUSY contributions in the most general way. In our analysis, 
we set the momentum transfer q 2 = 0. 

We first compute loop effects only (R-parity being conserved) by scanning over the MSSM 
parameters in the ranges shown in Table. HJ Here, tan/3 = v u jv& is the ratio of the up and 
down type Higgs vacuum expectation values; /i is bilinear Higgs coupling in the supersym- 
metric Lagrangian, which gives the mass to the Higgsino; Mi, M 2 and Mg are the masses for 
U(l)y, SU(2)l and SU(3) C gaugino, respectively; M\ are the diagonal mass parameters 
for the left- and right- handed squarks and sleptons of generation i, while M~ are the 

JLR 

left-right mixing parameters. In order to avoid unacceptably large flavor-changing neutral 
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Parameter 


Min Max 


tan (3 


1.4 60 


M 


50 GeV 1000 GeV 


M} 

Jt.r 


-1000 GeV 1000 GeV 



TABLE I: Ranges of SUSY parameters scanned. Here, M denotes any of \fj,\, Mi^,g, or the diagonal 
sfermion mass parameters M l ~ . The [i parameter can take either sign. The generation index % 

J L.R 

runs from 1 to 3. 



currents, we do not allow for flavor mixing between squark generations but do allow for 
superpartners of different generations to have different masses. 

In randomly choosing values for these parameters, we follow the conventional practice of 
using a linear distribution for all parameters except tan/3, for which we use a logarithmic 
distribution. We discard any points that yield SUSY particle masses below present collider 
lower bounds. In addition, we impose constraints from the Z-pole electroweak precision 
measurements, which are embodied in terms of the three oblique parameters S, T and U 
jl^j ]. To that end, we express the SUSY shifts in and k v in terms of these quantities: 



NC 
uN 



5 k,, 



a 



c 2 — s 2 I \ As 2 c 2 
crfWg 2 ) n Z7 (Mf 



S - aT + (J(t 



+- 
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+( 



VB 
SUSY 



c 2 — s 2 



II 77 (M|) Ad-iSvsY 



a. 



4g 2^SUSY^ 



(V 2 



(25) 



where remaining terms proportional to q 2 are negligible and have been dropped. Here A6 
is the SUSY contribution to the difference between the fine structure constant and the 

SUSY 

electromagnetic coupling renormalized at \i — Mz' Ad = \a{Mz) ' 



a 



Note that only S and T enter these expressions. Since these parameters are correlated, 
we use the 95% CL. S — T constraints |2j and retain only those parameter choices consistent 
with these constraints. We observe that this procedure is not entirely self-consistent, since 
we have not taken into account non-oblique corrections to Z-pole observables in deriving 
the oblique parameter constraints. Nevertheless, the essential, qualitative implications of 
precision Z-pole data for SUSY loop effects on the u^-q parameters are unaffected by this 
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inconsistency. We also omit parameter choices generating too large a SUSY contribution 
to the muon anomalous magnetic moment [1^]. Doing so limits left-right mixing for second 
generation sleptons to be fairly small. 

Before presenting our results, we also comment on the inputs used in Ref. jsj]. In that 
study, the authors only investigated the slepton contribution with four MSSM parameters: 
Mi, //, tan (3 and Mj, assuming the GUT relation for the gaugino masses and the slepton mass 
degeneracy. In our analysis, we include the contribution from both the squark and slepton 
sectors, taking into account sfermion left-right mixing, and allowing for non-universality 
between generations. 

With a sample of about 3000 randomly selected parameter sets, we calculated the MSSM 
contributions to R v and Rp. The numerical results are shown in Fig. Ufa). We observe that 
R v and R D are highly correlated. This correlation arises because MSSM contributions to 
the r^ l \g c ^ ) 2 term in Eq. (JIJ are small, while the contribution to the (gf 1 ) 2 terms are the 
same for both R v and R D . We also find that the magnitude of 8R VjD is dominated by the 
SUSY contributions to T (see Figs. EJb) andEJ, which is sensitive to isospin-breaking in the 
SUSY spectrum. It is bounded above by other precision electroweak data, thereby limiting 
the size of possible SUSY loop contributions to R v and R D to be considerably smaller than 
the deviations in Eq. (J2J). A breakdown of the various classes contributions is given in 
Fig. EJ More significantly, the sign of the SUSY loop corrections is nearly always positive, 
in contrast to the sign of the NuTeV anomaly. 

We do, however, find one corner of the parameter space which admits a negative loop 
contribution. This scenario involves gluino loops, whose effect can become large and negative 
when the first generation up-type squark and down-type squarks are nearly degenerate and 
left-right mixing is close to maximal. In this particular region of the MSSM parameter 
space, there is no gluino loop correction to the quark charged and neutral vector currents, 

2 3 

while the axial currents receive large corrections. When M~ = and and the second and 

JLR 

third generation soft parameters are sufficiently heavy (corresponding to a small value for 
T), the gluino contribution could give rise to a negative correction to R Uj p, as shown in 
Fig- El (a). In fact, the gluino contribution could be as much as fewxlO -3 in magnitude. 
As illustrated in Fig. EJb), however, equal and large left-right mixing for both up- and 
down-type squarks is inconsistent with the other precision electroweak inputs, such as the 
M\y and charged current universality . The shaded region in Fig. Efb) shows the region 



11 



x 10" 




x 10 




FIG. 4: Plot (a) gives the MSSM contribution to R u and R v with MSSM parameters chosen 
randomly from range shown in Table. [I] Plot (b) shows the dependence of 5R V on the oblique 
parameter T with the random MSSM parameter set. 
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FIG. 5: Contributions to 5R U from 5p NC — 5p cc (dashed line), 5k (dotted line) and (dash- 
dotted line). The solid line is the sum of all the contributions to SR V . Plot(b) shows the size of 
various components of 5p NC — 5p cc (solid line): T (dashed line), 5y B (dotted line) and 5yg — 5yg 
(dash-dotted line). The x-axis is the common first generation squark mass and first and second 
generation slepton mass. The other MSSM parameters are chosen to be tan (3 = 10, 2M\ = M2 = 
p = 200 GeV. The mass for the second and third generation squarks and third generation slepton 
are taken to be 1000 GeV. 
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FIG. 6: Plot (a) shows the MSSM contribution to R u (solid line) and R v (dashed line) with respect 
to the heavy first generation squark masses M~ = Mj = M~ , when the light first generation squark 
mass is fixed to be 200 GeV and the left-right mixing angle in fixed to be 9\ = 9 1 - = tt/4. No 
L — R mixing for other squarks and sleptons is assumed. The rest of the MSSM parameters are 
chosen to be: 2Mi = M 2 = n = 200 GeV, M- g = 200 GeV, Mj = 1000 GeV, M?' 3 = 1000 GeV 
and tan (3 = 10. Plot (b) shows the contours for the shifts 5R U (solid line) and 5Rp (dashed line) 
in the plane of slepton mass and heavy first generation squark masses. Slepton mass degeneracy 
is assumed. The rest of the MSSM parameters are the same as those used in plot (a). The shaded 
region in (b) shows the region preferred by charge current universality. 



preferred by these other inputs in the (M}) heavy and M~ t plane. Note that a shift in |V^ S |, as 
implied by the recent analysis of charged decays in the E865 experiment at Brookhaven 
National Laboratory [if]], would increase the shaded region in Fig. 

It is interesting to note that even if other precision data did not rule out this gluino loop 
effect, it could not account for the apparent deviation of sin 2 6 W from the SM prediction 
implied by the NuTeV analysis. The latter relies on a modified version of the Paschos- 
Wolfenstein relation 11 



R~ = 



R v - rR D 
1 — r 



5" 



2sin 2 6 W ) + 



(26) 



where we have shown the SM prediction for R with the + • • • indicating higher-order 
corrections. In practice, the extraction of sin 2 9 W relies on a slightly different combination 
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of R u and Bp in the numerator, 



Br = 



B v — £Rp 



(27) 



1 — r 

with £ chosen to be slightly different from r in order to minimize charm quark mass un- 
certainties in the sin 2 ^- extraction^]. It is straightforward to show that gluino loop 
contributions to B~ are proportional to £ — r, thereby minimizing their impact on sin 2 9 W . 
Specifically, we find that for maximal left-right mixing 9\ = 6^ = it /A, 

a s fJ . s 1 + r ^ /_ 5 



B~ 



— (Z-r) 



1 - -s' 
9 



2V 2 



M S , (Mi) heavy , (Mi) light 



V-2 



M~ g , (MJ) heavy , (Mi) heavy l - V 2 [M 5 , (Mi) %h t, (Mi) light 



(28) 



where a s is the strong coupling constant and V%{M, m\,m2) is defined in Eq. (IBlljl . 

In summary, it is difficult to explain the NuTeV anomaly in the i?-parity conserving 
MSSM framework. In most of the MSSM parameter space, the MSSM contributions to 
R v and R„ are small and have the wrong sign. Gluino loops could give sizable, negative 
contributions to B v and Bp, when the left-right mixing in the first generation squark sector 
is equal and close to maximum. However, this scenario, which is inconsistent with charged 
current data, could not in any case account for the deviation in sin 2 6 W implied by the NuTeV 
analysis. 



IV. RPV CONTRIBUTIONS 



As in the case of SUSY loop corrections, the effects of RPV contributions to B v and B 9 are 



12 1 . The relevant 



correlated with similar effects on other precision electroweak observables 
correlations are indicated in Table, m where we list the RPV contribution to four relevant 
precision observables: superallowed nuclear /3-decay that constrains \V u a\ 



measurements of the cesium weak charge j^], the e//i ratio R e /n in ratio 7T/2 decays 
and a comparison of the Fermi constant with the appropriate combination of a, Mz, 
2l|. The values of the experimental constraints on those quantities are given in 



19j], atomic PV 



and sin # p 



the last column. 

Given the experimental constraints on the first four quantities in Table m we obtain the 
one a allowed region for A 12 A;(e^) and /S! 2lk {d k R ) shown in Fig. [7| (a) by the solid ellipse. 3 



3 In performing the fit, we allow the signs of A^fe, A^. fe to be unconstrained. 
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A iifc(4) 

_L A.t\j \ IX' 


A ifci(<7£) A 




A 21fc(4) 


A 2fcl(4) 




s\v ud \ 2 /\v ud \ 2 


2 





-2 








-0.0029 ± 0.0014 


SQw/Qw 


-4.82 


5.41 


0.05 








-0.0040 ± 0.0066 


8R e /,j 


2 








-2 





-0.0042 ± 0.0033 










1 








0.00025 ± 0.001875 


5R U 








-0.21 


0.22 


0.08 


-0.0033 ± 0.0007 


5Rp 








0.077 


0.132 


0.32 


-0.0019 ±0.0016 



TABLE II: RPV contributions to 5\V ud \ 2 /\V ud \ 2 , SQ$/Q$, 5R e/fl , 5G^/G^, 5R U and 5R D . 
Columns give the coefficients of the various corrections A^- fe and A\2k to the different quantities. 
The last column gives the experimentally measured value of the corresponding quantity. 

Since the RPV corrections A-^ fc oc |A^ fc | 2 > 0, the physically allowed region - indicated by 
the shaded region in Fig. 0(a) - corresponds to all of the A\j k in Table ITD being non-negative. 
Taking the values of Ai2fc(e^) and A^^g^) from this region, we obtain the allowed shifts 
in R v and R„ shown in Fig. 0(b), dashed ellipse. We also show the corresponding 95% C.L. 
region (solid line in Fig. 0(b)). For simplicity, we have set A' 2kl (d l l) = since a non-zero 
value would yield only a positive contribution to these quantities. Even so, the possible 
effects on R u and R D are by and large positive. While small negative corrections are also 
possible, they are numerically too small to be interesting. 

We also performed a fit to the five RPV parameters including in addition the NuTeV 
results for R v and Rp. The one a allowed region for Ai 2 / C (e^) and A2 lfc (<4) i s given by the 
dashed ellipse in Fig. 0(a). At 95% C.L., at least one of the A-^ fc must be negative in this fit. 
Thus, inclusion of the NuTeV results appears to exclude the RPV SUSY effects summarized 
in Table ITTl with fairly high confidence. 

V. CONCLUSION 

The NuTeV anomaly remains in need of explanation. If one is ultimately unable to 
account for it with conventional, SM effects {e.g., small effects in parton distribution func- 
tions), then one would require an explanation lying outside the SM. Here, we have shown 
that such an explanation would be hard to come by in the MSSM alone. In general, SUSY 
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FIG. 7: Plot(a) shows the 1 a allowed region in A^^e^-A^^d-^) plane, with the best fit value 
denoted by the cross. The solid (dashed) ellipse is the fit excluding (including) the NuTeV R v ^ 
results. Shading indicates the physically allowed region, corresponding to A^ > and A- jfc > 0. 
Plot(b) shows the prediction for 5R U and 5Rp, using the 95% C.L. (solid line) or 1 a (dashed line) 
allowed values for Ai2fc(e^) and A 21fc (G^) from fitting to S\ V u d\ 2 /\V u d\ 2 , SQw /Qw > <^e/|* an d 
SG^/G^ with IS.' 2kl (d k L ) set to zero. 

loop corrections to R v and Rp generally have too small a magnitude and the wrong sign 
to account for the effect. An exception occurs for significant left-right mixing among first 
generation squarks, where a sizable effect of the necessary sign is generated by gluino loops. 
At present, precision data exclude this scenario with a high degree of confidence, though a 
new analysis of decays may weaken these restrictions considerably. Even in this case, 
however, the value for the weak mixing angle extracted from neutrino-nucleus scattering will 
be largely unaffected by gluino loops when a Paschos-Wolfenstein type relation is used for 
the extraction. 

R parity-violating effects also appear at odds with the NuTeV anomaly. Inclusion of 
these effects could resolve an apparent conflict between tests of charged current universality 
and implications of SUSY-breaking models for the MSSM spectrum 9|]. However they would 
also render the LSP unstable and, therefore, rule out SUSY dark matter. At face value, the 
NuTeV results appear to disfavor this resolution of the charged current universality problem. 
In short, we conclude that the MSSM - with or without R-parity conservation - is likely 
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not responsible for the NuTeV anomaly. The culprit, apparently, is to be found elsewhere. 
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APPENDIX A: ESSENTIAL FEYNMAN RULES 



The complete set of Feynman rules for the MSSM can be found in |13|, |22|. Here, we give 
only a brief list of relevant vertices. 

The fermion-sfermion-gaugino vertices are: 

Fr , Fj , Q P 



Fl 1 Ft 1 Qf 

Xp or Xt c * X° P ~g a 

9l } ' P Pl + 9 F / P P R ) ie [itP L + g!tP R ) ig sKp (0P L + g^P R ) 

where e is the absolute value of the electron charge, gs is the SU(3) C coupling constant, and 



X^p are Gell-Mann matrices j22|, normalized according to tr[X a X b ] = l/25 ab . We use the 
capitalized letters I and J to denote the family index for quarks and leptons (I, J = 1, ■ ■ • , 3), 
small letters i and j to denote the index for squarks and sleptons = 1, • • • , 6 except for 
sneutrino, when i,j = l,---,3), and small letters p and n to denote the index for the 
neutralinos (p,n = l,---,4) and charginos (p,n = 1,2). The index a is reserved for the 
gluino index, a = 1, • • • , 8. 

The first vertex represents the coupling of the fermion Fj to the sfermion F[ and chargino 
Xp (or its charge conjugation field Xp C it F = D, L.). The chargino coupling constants g^f^ 
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for the quark sector follows (the repeated index is summed over): 



U Tip 

9l 



Ujip 

9r 



D jip 

9l 



Djip 

9r 



-V, 



i j 



-Vr 



z%u; x 



muj 



m Dl 



u~ 



•s " V2M W sin /3 

7*1 



S 



V2M W cos p 

ZpUp2 



U p2 



-Vi 



IJ- 



nr*lllf* ■■ -u I y*I+3, ii r* 

Zu Vpl ~ V2M w smp Zu Vp2 y 
m Dl 



(Al) 



s V2M w cosp^ U p2 ' 
where s (c) is the sine (cosine) of Weinberg angle, Vu is the usual CKM matrix, m^m^ 
is the mass of "up" or "down" quark for generation I, tan/3 = v u /vd is the ratio of the 
expectation values of the Higgs scalars, Zjj t u are the unitary 6x6 mixing matrices for up 
and down squarks, respectively, that diagonalize full sfermion mass matrices, and U and V 
are the unitary mixing matrices that diagonalize the chargino mass matrix [l3| . In practical 
calculations masses of the first and second generation quarks can be neglected as they are 
much smaller than the mass of the W boson. Similarly, for lepton sector, 



V Tip 

9l 



9r 

Ljip 

9l 



9r JW 



-A z lu;i 



m Lj 



V2M W cos (3 



ZJ +3,2 j r* 
r. U, 



p2 



v V pl 
S F 

m Lj 



Zt J %2, 



(A2) 



s y/2M w cos[3' 

Note that the mixing matrix for sneutrinos Z v is 3x3 because in the MSSM neutrinos are 
purely left-handed. 

The second vertex represents coupling of a fermion to a sfermion and a neutralino. The 
corresponding coupling constants g [ ip R for the fermion-sfermion-neutralino vertex are: 



Ujip 

9ol 



Ujip 
9 R 



Djip 

9ql 



Djip 

9or 



9ol 



1 



\/2sc 
1 /4 



Z*Ji 
U 



n; 2 c + -n; iS ) - 



m Uj y *J+3,i AT * 

M z sinp Zu 7 > 4 



muj 



_ y -l ■ '■>. ': \ _ 

V2c\3 U pl M z smP 



Z v Ji N pi 



\/2sc 
1 



ZJi 
D 



M z cos (3 



1 Z£ +3 > 1 N, 



P 3 



2 7 J+3,i AT . m Dj 
-Zj D iVpi + — 

3 M z cos p 



y/2c 



V2sc 



p2' 



n; iS 
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9m 

Ljip 

9ol 



Ljip 

9m 







V2sc 
1 



zi 1 (n; 2 c + n; iS 



m Lj 



7j+3,i 



M z cos (3 



N* 



P 3 



2Z 



J+3,i 



N pl + 



:Z J L l N p3 , 



V2c V u " P± ' M z cos /T 
where N is a 4x4 mixing matrix that diagonalizes neutralino mass matrix 
Finally, quark-squark-gluino (the third vertex) couplings 9gl gr are: 



(A3) 



tit- 

Dji 

9gl 



-y/2Zg 



9g j r = V2Z, 



*J+3,i 
(7 



(A4) 



The gauge boson-gaugino-gaugino couplings are given as follows: 




W-" 




Z'-' 



^7" (0%Pl + 0%Pn) i%f (0%P L + 0^'P R 




where 

o L pn = (n p2 v: x - ^N pA v:)j eft = [N; 2 u nl + -^iv; 3 c/ n2 J 

£' ( t/ t/* ~^t/ t/* i £ jA nH' /" rr* rr ^ rr* rr i x „2 



R = 

^pn 


[N; 2 u nl + 


o A = 

^pn 


{-u; x u nl - 


R " = 

^ pn 


^pn • 



o L p : = (-\n p3 n; i3 + \n pA n:^ o% = -o;t". (as) 

In addition, for the photon-chargino-chargino vertex A^ 1 — xt ~ X p is given by —ie , y IJ, 5 pn . 
The gauge boson-sfermion-sfermion vertices are given by: 



F / F,; 



PF t/ P Fi 



/ 



W ±fi 'wvw V = Z' 1 

J \ ~ J \ 

v F' \ p. 

ie 9 W m + P^)" " /A/./ •/'/; + PR-)'' 
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For the W- coupling: 

9wm = -^ s VjlZ ^ Z u\ QwUv, = -^Z^ZfK (A6) 
and for the Z and 7 couplings: 

9zf z Fj — ~— (l[Z F Z F Ij - Q F s 2 5ij^ 

g lhFj = -Q F S ij . (A7) 



Here, 7-f and Q F are the isospin and the electric charge of the sfermions F, respectively. 

The two scalar-two gauge boson vertices necessary for the calculation of the gauge boson 
self energies are (see Ref. |22j]): 

F F: 

v y -2 



V' 1 vwvwvWv y' v 



1 1 j 



where 



n - - — 7 Ki 7 *Kj ny 

yWWFiFj — 2 S 2 

g zzhh = ^- 2 {2li{li -2Q F s 2 )Z^Zf> + 2Q F s%)g 

J s c x ' 

g ZlhFj = -Qf (2l[Z^Zf> - 2Q F s%) gT 



sc. 

= 2Q 2 F 5ijg^ u . (A8) 



APPENDIX B: EXPRESSIONS FOR THE RELEVANT FEYNMAN DIAGRAMS 

In this section, we list analytical expressions for the MSSM contribution to the gauge 
boson self-energies, the fermion wave function renormalization, the gauge boson-fermion- 
fermion vertex correction, and the box diagrams relevant to the neutrino-nucleus scatterings. 



Modified dimensional reduction (DR) scheme is used in computing MSSM loop contribu- 
tions, although in the Appendices, we have neglected the hat in all relevant variables. 

1. Gauge Boson Self-Energies 

We first define the following class of two-point integration functions: 

,k„) = />MH + ( 1 - i k- I (.- i »44 m 
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ft 



(b) 



V vwvwVwv^' 
(c) 



FIG. 8: Feynman diagrams of one loop SUSY contributions to the gauge boson self energy. 

F 12 (m 1 ,m 2 ,m 3 ) = F 1 (m 1 ,m 2 ,m 3 ) - F 2 (m 1 ,m 2 ,m 3 ), (B2) 

where /i is the renormalization scale. In our analysis, m\ will be replaced with either the 
external fermion mass squared or one of the Mandelstaam variables. For the first two 
generation quarks and the first two generation leptons the fermion mass can be set to zero. 
Focusing on the case when all Mandelstaam variables are small compared to m\ 2 , we obtain: 



i 7 o(mi,m2,0) = him\ — 1 



rrin 



rrin — m 



In — 2 — In /i 



F 1 (m 1 ,m 2 ,0) 



2\nm\ - 1 + 



2ml 



2m\ 



— 1 In — - 

m\ — m 2 (m 2 — mf ) 2 m\ 



- 2 In /i 2 



(B3) 



a. Gaugino Loops 



The Feynman diagrams of this type are shown in Fig. |Sfa). The contribution of gaugino 
loops to W and Z self-energies is given by 

^vvil 2 ) = "f E W + 9L9R)m Xn rn Xp F (m Xn ,m Xp ,q) 

^ p,n { 

+ (\9l\ 2 + \9r\ 2 ) [2q 2 F 12 (m XnJ m Xp ,q) - m Xn F x {m Xn ,m Xp ,q) - m^F^m^, m Xn , g)] |(B4) 

where g M is the momentum carried by gauge boson V, and q 2 = q^q a - The couplings gi and 
gR are listed below for Tlww and Hzz- 





Xp 


Xn 


9l 


9r 


comment 


n x°x + 
LL WW 

TTX+X+ 

LL zz 

n x°x° 

lL ZZ 


Y° 

A.p 

Xp 


Ap 


Xn 

Xn 


An 


O l p Js 

L Jsc 

O^/sc 


0™/sc 


multiply 1/2 
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The Z — 7 mixing tensor and photon self-energy from chargino loops is given by 



7T p 

with the couplings gy and gy listed below: 



(B5) 



Uvv 


9v 


0v/ 


comment 


TT X+ 

77 


(oi; + o%)/sc 
-i 


-1 
-1 


multiply 2 



6. Scalar Loops 

Diagrams of this type are shown in Figs. EJb,c). The contribution of Fig. |H] (b) to Ilvi/w, 
Tlzz, n^ 7 and II 77 is given by 

-,2 



n^(g 2 ) = ~Y,9VF i F i 9v'F i F j \ m % + 



2 Q 

m f 



nr 



'?.Fi(mp ,mp q) +m%F 1 (mp, ,mp ,q) - q 2 F l2 (nip m p q) 



,(B6) 



The sfermion pair (Fi, Fj) running in the loop for Uvv is listed below: 



(V, V) 


(FuFj) 


(W,W) 


(h 


(Di,Uj) 


(Z,Z) 




Lj), (U h Uj) ,(Di,Dj) 


{Z,l) 


(Li, Lj) 


(Ui,Uj) ,{D h Dj) 


(7,7) 


(Li, Lj) 


(Ui,Uj) ,(A,A) 



For squark contributions, an additional color factor of N c enters the right-hand side of 
Eq. flBHJ). 

The contribution of Fig. |H] (c) to the self-energies is given by 



u vv(<l 2 ) = ^J29vv>FF m % (l - lnm % + ln K 



(B7) 



where F = v, L,U and D except for V = 7, for which v does not contribute. 
When q 2 — > 0, Uz-y(q 2 ) / q 2 and Ii 11 (q 2 )/q 2 from sfermion loops reduce to 

rw(<? 2 ) 



a 



Y,9 V F i F i 9v'F i F i ( ln77l % ~ 1] V) ■ 



(B8) 
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FIG. 9: Feynman diagrams of one loop SUSY contribution to the fermion wave function renormal- 
ization. 

2. Field Strength Renormalization for Fermions 

The diagrams contributing to the field strength renormalization for fermions all have the 
form shown in Fig Therefore, the same formula can be applied to contributions with 
charginos and neutralinos in the loop, provided that the appropriate couplings and masses 
are used. The gluino contributions are multiplied by an extra factor of Casimir factor 
C*2(iV) = 4/3 for SU(3) C . The field strength renormalization for the left-handed quark Fj is 

i,p i,p 

+t x ^J2\9GL\ 2 Fi{mp.,mg,m Fj ), (B9) 

i 

where m- g is the gluino mass and F' stands for the isopartner of the fermion F {e.g. if F is 
the up-type quark, then F' is the down-type quark). By replacing L — > R one easily obtains 
the field strength renormalization for the right-handed fermions. The same formulae apply 
to the sleptons provided that the appropriate couplings and masses are used. Naturally, 
corrections involving the strong coupling are absent in that case. 



3. Vertex Corrections 

One-loop SUSY corrections to the V — f — f vertex are shown in Fig. There are 
two types of corrections: loops with the vector boson coupling to the scalar particles, Fig. 
HUT a), and loops with the vector bosons coupling to the gauginos, Fig. ITuT b). The complete 
set of Feynman diagrams for each vertex, such as W — d — u, can contain more than one 
diagram of each type. However, since these diagrams differ only in the specific values of 
the masses and couplings needed to obtain the numerical answer we do not show all graphs 
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FIG. 10: Feynman diagrams of one loop SUSY contribution to the gauge coupling vertex. 



explicitly. Below, we use the superscripts (a) and (b) to distinguish between the two types 
of loop graphs. To distinguish the graphs of the same type, we supplement the superscript 
with a number. 

We first define the three-point integration functions needed for the evaluation of the 
triangle diagrams: 



V\{M, mi,m 2 ) 

V 2 {M, mi, m 2 ) 
D 3 (M,m 1 ,m 2 ) 



I dx dy 
o Jo 



y 



D 3 (M,m 1 ,m 2 ) 

1 f 1 r 

dx / dyylu D 3 (M,m 1 ,m 2 )/(i' 
o Jo 1 

(1 - y)M 2 + y[(l - x)m\ + xraj]. 



(BIO) 



Explicitly: 



Vi(M } m l} m 2 ) 
V 2 (M,m l} m 2 ) 



m 2 ln^ 



(M 2 — m 2 )(m 2 — m 2 ) (M 2 — m 2 )(m 2 — m 2 ) 
1 



21nM 2 - 3 + 



2m* 



ln m? 



(M 2 — m 2 )(m 2 — m 2 ) M 2 



2m^ 



In 



m 



{M 2 - m 2 2 ){m{ - m 2 2 ) M 2 
Defining the matrix element for the vertex to be 



2 In// 



(Bll) 



(B12) 



Diagram (a) gives 



rL;(a) 
VF T F, 



If J 



rR;(a) 
VF T Fj 



a 
a 



E 9 VFiF, 9xl9 J xL V 2 {™ X p > m P, , m i 



2VJ 



~TzJ2 9vF 1 F j 9 x R9xR V 2( m x^ m F^ m i 

P,i,j 



(B13) 
(B14) 
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Diagram (b) gives 



°VF I F J 



,R;(b) 
'VFzFj 



a 
Air 
a 
Air 



p,n,i 

£[■ 

p,n,i 



'XX^xLyxL '"-Xp'^Xn^l yVxx^xLyxL 

9vxx9 X R9xR m x P m Xn V i ~ 9v xx 9 x r9 x r (1/2 + V 2 ) 



(B15) 
(B16) 



with the argument for functions Vs in Eqs. f|B15|) . ()B16|) being V 1)2 (m^, m Xp , m Xn ). The 
explicit expression for gyp.p., 9v xx i 9 x l an d 9 x r for each individual vertex diagrams will be 
given below. 



a. Charged Current Vertex 
The vertex correction to W — Dj — U,j is 



°VF,Fj 


Xp -Pi -Pj' 


9vF i F j 


9 x l,R 


n J 

9 X L,R 


S (b) 

"VFjFj 


Xp Xn Fi 


9v X x 


9v X x 


9 X L,R 


9 X L,R 


x(«l) 


X° A 


9wD i U j 


Djip 
9oL,R 


Ujjp 
9oL,R 


r(61) 
"WDjUj 

°WD I Uj 


Xp Xn -A 
Xp Xn Ui 


Otis 


-o:i/s 


Djip 
9oL,R 

Diip 

9l,r 


Ujin 

9l,r 

JJjin 
9oL,R 



The vertex correction due to gluino exchange is similar to S\^ DlU , with the substitution 



of 



a 



as, 



Xp 



~ Djip 

9-, 9ol.r, 



Dji 

9gl,ri 



Ujjp Ujj 

9ol,r ~^ 9gl,r, 



(B17) 



and multiplication of the whole expression by 4/3. This substitution rule also applies for 
the neutral current vertex listed in the next section. 

Similar to the case of the field strengths, the corresponding vertex corrections involving 
leptons SyLjuj are obtained from by using the appropriate masses and couplings. The gluino 
loop corrections must be omitted in that case. 



b. Neutral Current Vertex 



For the Z — Dj — Dj vertex 



°VF I F J 


Xp Fi Fj 


9vF i F J 


9 x l,R 


9xL,R 


S (b) 

"VFjFj 
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9v X x 
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°ZD I D J 
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°ZD I D J 
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Xt Ui Uj 


9zbibj 
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Djip 
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9l,r 


Djjp 
9oL,R 

Djjp 

9l,r 


r(M) 
°ZD I D J 
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°ZD I D J 


Y° v° D 
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-0 L p Usc 
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9l!r 



25 



For the Z — Uj — Uj vertex 



x(a) 


X P ^ Fj 




9 X L,R 


n J 


°VF I F J 


Xp Xn F{ 


9v X x 


9v X x 


9xL,R 


9 X L,R 


r(al) 
°ZV I V J 
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xl Ui u 3 
xt A A 


9zUiU 3 


Ujip 
9oL,R 

Uiip 

9l,r 


Ujjp 
9oL,R 

Ujjp 

9l,r 


r(W) 
°ZU I U J 

°ZU I Uj 


xl xl Ui 

Xp Xn A 


0%/sc 
O^p/sc 


0%/sc 
0%jsc 


Ujip 
9oL,R 

Uiip 

9l,r 


Ujin 
9oL,R 

U jin 

9l,r 



The radiative corrections to the lepton neutral current vertex are directly obtained from 
the above expressions by replacing the up (s)quark with the (s)neutrino and the down 
(s)quark with the (s)electron. 



4. Anapole Moment Corrections 

In the presence of parity-violating interactions, higher-order contributions can generate 
the photon-fermion-fermion coupling of the form (see, e.g. Ref. 

<// = - e ^f^V~ ^)7 5 /^- (B18) 
The contributions from gaugino-sfermion loop at q 2 — > are: 
FaM) = F$tf )+F W {q * )+F V {q i) 



Ft/if) = -Qf^B y £[\9S\ t -\9& 



487r v " u ^ 1 J Jo (1 - x)m 2 + xm 2 f 

l iP Xp Ji 

Oil\/L y ^ — -v /. f Vri 1 9 i fi'ni r >\ X^ (Ix 



48/T f-f 1 l3i< 1 ) Jo (I — x)m 2 + + xm\ 

l 'P K ' fi Xp 

where /' stands for the isopartner of the fermion / {e.g. if / is the neutrino, then /' is the 
electron). Notice that for quark anapole moment, an additional gluino contribution should 
be added, when parity is broken in the presence of a non-zero left-right mixing in the squark 
sector and non-equal diagonal left and right squark masses. The gluino contribution can be 
obtained from F^f using the substitution rule give in Eq. (|B17|) . It should be remembered 
that the anapole moment of an elementary particle is not a physical observable j^. Here, 
we separate it from the other one-loop contributions purely for clarity of presentation. 
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5. Box Graphs 

Let us introduce the following notation: 

Hf'J)x = /V)7a(1-75)/(p) 

R(f, f)x = /V)7a(1 + l5 )f(p), (B20) 

where p and p' are the momenta of the incoming particle / and outgoing /', respectively. 
The four-point integration functions are defined as: 

B l (M 1 ,M 2 ,m 1 ,m 2 ) = C dx f dy f dz n ,~ Z * r 

Jo Jo Jo -L> 4 (Mi, M 2 , mi, m 2 ) 

rl rl rl zCl — z) 

B 2 {M 1 ,M 2 ,mi,m 2 ) = / dx dy dz 

Jo Jo Jo D%{M 1 ,M 2 ,m 1 ,m 2 ) 

D 4 (M 1} M 2 ,m u m 2 ) = z[(l - x)M\ + xM\\ + (1 - z)[(ym\ + (1 - y)m 2 2 \. (B21) 
The explicit formulae for Bi and B 2 are 

m\ ln^| 

B\{Mi, M 2 ,m 1 ,m 2 ) 



2{m\ - Ml){m\ - m 2 2 ){ml - Mf) 
mfln^f Mflnf 



+ 



2(m| - m\){m\ - Ml){m\ - Mf) 2(M X 2 - m 2 )(M 2 - m 2 2 ){M\ - Mf) 

2 i A?"? 

mf in — §■ 

B 2 {M 1: M 2 , mi ,m 2 ) 



(m{ - Ml) {ml - m|.)(m 2 - M|) 
ml In ^ M^ln^ 



' (m| - m?)(m| - Ml) (m| - M|) + (Mf - mf ) (Mj 2 - m|)(M x 2 - M|) " (B22) 

a. Charged Current Box 

We define the matrix element of the box diagram as 

Mg c x = -&5° b c L{im, u,) x x L(u, d) x + ... (B23) 

to factor out the overall Fermi-constant. The dots denote other Dirac structures appearing 
in the amplitude that make negligible contributions (suppressed by m^) to the total cross- 
section. 

CC;a . cCC;b . x CC;c . x CC;d 



For the graphs in Fig. [TT]we have: 5 B C = 5 B ' a + 5 B ' + 5 B ,c + 5 



B 



S b = 1^ 9ol 9l 9ol9l m xl m xi B ^ m ^ m D^ m x^ m xi' 
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FIG. 11: Feynman diagrams of one loop SUSY contribution to neutrino charge current box diagram. 
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Similarly, the box diagram contributing to the muon decay \x — > Vpev e can be obtained from 
5% using the substitution: 



5% c — > 5^*, d — > e, u — * i> e , — ► Zj, £7^ — > 9j, 

djp ejp ujn u e jn 

9l,r,ol,or ~^ 9l,r,ol,or-> 9l,r,ol,or ~ * 9l,r,ol,or- 
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b. Neutral Current Box 



We define the matrix element of the neutral current box diagram as 

Ml C x = [S&Lfa v,) x x L(q, q) x + 5§ q L(fi, v,) x x R(q, q) x 

The explicit expressions for the Vu — q box graphs are given below. 
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Up-quark box diagrams (Fig.ITJ: = 5^ u ' a) +S^ u ' b) + 5^ c \ 5% u = +5? ( * h) 
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FIG. 12: Feynman diagrams of one loop SUSY contribution to v„ — u neutrino neutral current box 
diagram. 
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FIG. 13: Feynman diagrams of one loop SUSY contribution to — d neutrino neutral current box 
diagram. 
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Down-quark box diagrams (Fig. fl~2|) : 5^' d = S^ d ' a ^ + 5^ + 8^ , 5 B ' d = S^ d ' a ^ + 

rR;(d,b) 
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APPENDIX C: RADIATIVE CORRECTION TO NEUTRINO-NUCLEUS 

INTERACTIONS 

Given the expressions for the fermion field strength renormalization SZ Rl R , vertex correc- 
tion 5yp [Fj and box diagrams 8g R , 5y R in Eq. (JHJ), S v in Eq. (JO]) and Sy R,q in Eq. (fTT|) can 

be expressed as 

Kb = \ {Ml + 6Z2 + 5Zt + 6Z?) - V2s (6^ e + 5^) + 8% (CI) 
5vb = ~ + Mi + 5ZI + 5Zfj - V2s (b L Wv ^ + 5^ du ) + 5g c (C2) 
5 V = 5Z? - 2sc 5 L Zv ^ (C3) 



tvT u ' d = \SZI - sc 6 L Zqq + 5^ (C4) 
*?T M = ~^ - sc 6 R qq + 6 R «. (C5) 

The neutrino anapole moment contribution to the neutral current neutrino-nucleus scat- 
tering is shown in Fig. [IJb), which is absorbed into 5 

5k = -4c 2 F A , u (q 2 ). (C6) 
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